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ABSTRACT 

We have performed power spectral analysis of surface temperatures, velocities, and magnetic fields, 
using spectro-polarimetric data taken with the Hinodc Solar Optical Telescope. When we make power 
spectra in a field-of-view covering the super-granular scale, kinetic and thermal power spectra have 
a prominent peak at the granular scale while the magnetic power spectra have a broadly distributed 
power over various spatial scales with weak peaks at both the granular and supergranular scales. To 
study the power spectra separately in internetwork and network regions, power spectra are derived 
in small sub-regions extracted from the field-of-view. We examine slopes of the power spectra using 
power-law indices, and compare them with the unsigned magnetic flux density averaged in the sub- 
regions. The thermal and kinetic spectra are steeper than the magnetic ones at the sub-granular scale 
in the internetwork regions, and the power-law indices differ by about 2. The power-law indices of the 
magnetic power spectra are close to or smaller than -1 at that scale, which suggests the total magnetic 
energy mainly comes from either the granular scale magnetic structures or both the granular scale 
and smaller ones contributing evenly. The slopes of the thermal and kinetic power spectra become 
less steep with increasing unsigned flux density in the network regions. The power-law indices of all 
the thermal, kinetic, and magnetic power spectra become similar when the unsigned flux density is 
larger than 200 Mx cm~ 2 . 

Keywords: Sun: atmosphere - Sun: photosphere - Sun: granulation - Sun: surface magnetism 



1. INTRODUCTION 

On the solar surface, interaction between magnetic 
fields and surface convection produces varieties of struc- 
tures over the broad spatial scale from 10 5 km, such 
as sunspots and active regions, down to ~ 100 km or 
even below, such as small flux concentrations. Spatial 
power spectra of velocity and magnetic fields on the so- 
lar surface provide a clue for understanding at which 
scale kinetic and magnetic energies are injected, trans- 
ferred, and dissipated on the solar s urface (see recent 
revie ws by iRieutord fc Rinconl 120101 : lAbramenko et al.l 
120121) . Basic information of the granular convection was 
retrieved from intensity and Doppler-shift images show- 
ing temperature and velocity fluctuations on the solar 
surface. Their power spectra were analyzed by various 
authors using high reso lution observations taken with 
ground-based telescopes llFrenkiel fc Schwarzschildlll955l: 
Beckers fc Parnelll Il96<l iReilingj 119711: IDurr ant fe Nesis^ 



1982 1 IMullerl ll989HRoudier et al.l 119911: lEspagnet et al.l 
1993T) . They discussed how the turbulent convection on 



the solar surface makes slopes of the power spectra to- 
ward a high wavenumber domain by comparing them 
with the Kolmogorov's fc~ 5 / 3 power-law, where k is the 
wavenumber. Power spectra of surface magnetic fields 
were also analyzed using longitudinal magnetograms 
taken with ground- b ased t elescopes and SOHO /MDI 
(iNakaeawa fc Priestj 119731: iNakaeawa fc Levind 1 1974 : 
Knobloch fc Rosnerll98lULee et al IU997UPetrovavl2001 



lar and sub-granular scales because of the lack of spatial 
resolution in polarimetric observations. 

In most studies using ground-based telescopes, power 
spectra were significantly affected not only by the in- 
strumental resolving performance but also degradation 
of the image quality due to atmospheric seeing, espe- 
cially when we study the power spectra at the scale 
smaller than granules thoug h calibration of such degra- 
dation was attempted (e.g. iLee et all Il997f ). It is also 
important to distinguish between different observed re- 
gions when we make power spectra, because unsigned 
magnetic flux densities have significant variation depend- 
ing on internetwork, network, and active regions. Influ- 
ence of magnetic fields to the surface convection can- 
not be negligible in strong field r egions while i t may 
be small in weak field regions. | Abramenkol (2005) 
and lAbramenko fc Yurchvshvnl (|2010l) studied magnetic 
power spectra in various active regions, and suggested 
correlation between the power-law index of the magnetic 
spectrum and flare productivity of an active region. But 
it has not been understood yet how the power spectra 
changes from internetwork regions to net work ones. 

Th e Solar Optical T elesco pe (SOT iTsuneta et all 
20081: iSuematsu et all 120081: llchimoto et al.l I200S 



Shimizu et al.ll2008D aboard Hinode ()Kosugi et al.ll2007 



Abramenko et all 120011 ). But it has been hard to study 
the power spectra with enough confidence at the granu- 

3 Dept. Astrofi'sica, Universidad de La Laguna, E-38206 La 
Laguna, Tenerife, Spain 



is a suitable instrument for studying the behavior of 
the power spectra at the granular scale and even below 
owing to the diffraction-limited perfo rmance of the 
50 cm aperture and stable image quality. IRieutord et all 
(2010) examined power spectra of intensities, transverse 
and longitudinal velocities using filtergram data taken 
with the Hinode SOT. We use two-dimensional maps 
taken with the spectro-polarimeter (SP) of SOT, which 
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not only provides intensity and velocity maps, but 
also magnetic field maps with the same instrument. 
It is becoming important to understand properties 
of small-scale magnetic structures in internetwork re- 
gions because the Hinode SP observations suggested that 
small-scale fields are possibly genera ted and sustained by 
turbul e nt small-scale dynamo (e . g. Ilshikawa fe Tsuneta 
2009at iPietarila Graham et al.l l2009t iDanilovic et al. 
20101: ILitesll201in . The most common indicator is the 
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probability distribution functions (PDFs) o f mag netic 
field strength or magn etic flux. iStenflol (|2012j ) and 
lAbramenko et all (|2012f ) used the magnetic power spec- 
tra for investigating the small-scale magnetic structures 
using the Hinode SP data. 

In this article, we make power spectrum analyses of 
the surface temperatures, velocities, and magnetic fields 
using the Hinode SP data, including calibration of the in- 
strumental resolution performance. Section 2 describes 
how the SP data are processed to make the power spec- 
tra. Section 3 presents the power spectra covering both 
the supergranular and granular scales. In Section 4, the 
power spectra are created in small regions to distinguish 
internetwork and network regions. We study how the 
power spectra change with the unsigned magnetic flux 
density. 

2. DATA ANALYSIS 
2.1. Data sets and reduction 

We use a data set taken with the SP of the Hinode 
SOT. The data set consists of 30 SP scans taken between 
November 2006 and December 2007. The SP recorded 
full Stokes spectra of the two Fe I lines at 630.15 nm 
and 630.25 nm with the normal map mode whose spatial 
sampling was about 0.15" per pixel and integration time 
is 4.8 sec per slit position. The field-of-views (FOVs) of 
the scans were equal to or more than 1000 steps in the 
scanning direction and 1024 pixels along the slit, and con- 
taining the solar disk center. It took about 85 minutes 
to complete the scan of each FOV. Their observing tar- 
gets were mostly in quiet regions, but some of the scans 
included active regions also. We identify 23 scans out 
of 30 containing only the quiet regions using the condi- 
tion in which an averaged density of unsigned flux {\Bz\) 
in the FOV is smaller than 20 Mx cm" 2 . The SP data 
are calibrated with the standard routine spjprep available 
under the Solar Software (SSW). The wavelength posi- 
tions of the two Fe I lines are calibrated with an average 
line-center position outside of a sunspot. 

The SP data provide two-dimensional maps of surface 
temperatures, line-of- sight (LOS) velocities, and mag- 
netic field vectors. The temperatures T are obtained us- 
ing continuum intensities averaged in the spectral win- 
dow outside of the Fe I lines and under assumption of 
black body radiations. We assume that average temper- 
ature T is 6520 K outside of sunspots to convert the 
observed continuum intensities to the temperature. The 
LOS velocities v z are obtained using Doppler shifts of 
the Fe I line at 630.15 nm. The magnetic field vectors 
(B x , By, B z ) are determined from wavelength-integrated 
Zeeman-i nduced pola r izatio n signals with the method de- 
scribed in lLites et al.l (|2008l ). The x, y, and z components 
correspond to ones in the scanning, slit, and line-of-sight 
directions, respectively. It is not necessary to convert 



-§ 60 



200 



400 




(Mx cm-2) 
100 



100 



40 60 
(Mm) 

Figure 1. Two dimensional maps of (a) the temperature fluctua- 
tion AT and (b) the vertical component of magnetic flux densities 
B z obtained with the Hinode SP on 24 April 2007, as an example 
of the data set used in the power spectrum analysis. The FOV 
is 1000 steps in the scanning direction and 1024 pixels along the 
slit, which covers a square of about 110 Mm. The hours of the SP 
scanning observation are shown at the top of the figure. 



the observed vectors into the ones in the local reference 
frame because the maps were taken near the disk cen- 
ter. To derive the transverse component (B x ,B y ), the 
180 degrees ambiguity in the azimuth angle has to be 
resolved. This is simply done by choosing the azimuth 
to be closest to the computed potential field using B z 
as an input. Fig. [T] shows two dimensional maps of the 
temperature fluctuation AT = T — Tq and the vertical 
component of magnetic flux densities B z as an example 
of the data set. 

2.2. Derivation of power spectra 

For deriving spatial power spectra, we apply the two- 
dimensional Fourier-transform to each map of the above 
parameters, 



P{f'i k x , ky) 



jp /( x ' y) ex P {-IniikxX ■ 



k y y)) 



x,y 



where f(x,y) is either AT = T — To, v z , B x , B y , or 
B Zl and N the number of pixels in one direction and 
1024 in our analysis. The wavenumbers in the x and y 
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directions are represented as k x and fcj,0 To have 1024 
pixels in both the x(scanning)- and y (slit)- directions for 
simplifying the Fourier-transform in Eq. (1), we extract 
1024 steps closest to the disk center when the number of 
steps in a scan is more than 1024. Because our selection 
criteria of the SP data set is that the number of pixels in 
the x-direction is more than 1000, the number of pixels 
in the x-direction is fewer than 1024 pixels in some of 
the maps. For those maps, we apply padding at both 
the ends in the x-direction to make a map consisting 
of 1024 x 1024 pixels. We confirmed that the influence 
of the edges and padding is negligible when we apply a 
Hamming window in the Fourier-transform. 

The two dimensional spatial power spectra P(f; k x , k y ) 
is then converted to their one dimensional form P(f; k) 
by integrating the two-dimensional power spectrum in 
azimuthal angle of the wavenumber (k x , k y ), 

k) = 2?rfc (P(f; k x , fci/)/fc'e[fe-Afc/2,fc+Afe/2) > ( 2 ) 



where k! = w k x + ky, and Afc is the sampling interval 

of the wavenumber, Ak = l/(NAx), where Ax is the 
pixel scale and 0.15" or 0.11 Mm on the solar surface. 
The power spectra thus derived are normalized using the 
following equations to represent them in a unit of energy 
density per unit wavenumber, erg cm -3 (1/Mm) , 



E th {k): 
E vz {k)-- 
E bh (k)- 
E bz {k)- 



3 nk B P{AT; k) 



2 T 



Ak 



1 P(v z ;k) 
2 P Ak 

1 (P(B x ;k) + P(B y ;k)) 
8tt Ak 
1 P(B z ;k) 



8tt Afc 



(3) 
(4) 
(5) 
(6) 



respectively, where ks is the Boltzman constant, n the 
number density on the solar surface, where we assume 
n = 1 X 10 17 cm~ 3 as a typical photospheric number 
density, p the mass density on the solar surface, where 
we assume p — 2.7 x 10~ 7 g cm~ 3 . The method used 
here is essentially c onsistent with the on es used in previ- 
ous studies (such as lRieutord et al.ll2010D except the nor- 
malization. The power spectra E t h(k), E vz (k), Ef,h(k), 
and Eb z (k) are created from each map. The power spec- 
tra averaged in the quiet regions (i. e. using 23 maps 
out of 30) are shown as dashed curves in Fig. [2] 

Because the power in the higher wavenumber domain 
is affected by the instrumental spatial resolution, it is 
necessary to calibrate the optical transfer function to ob- 
tain the shape of the power spectra on the solar surface. 
The modulation transfer function (MTF) of the instru- 
ment is evaluated at 630 nm and is shown in Fig. [3l 
Though the telescope achieved diffraction limited perfor - 
mance (|Suematsu et alJl2008t iWedemever-Bohml l2008h . 

2 The definition of the wavenumber is the same as the one in 
Ricutord ct al. (2010), in which the spatial scale corresponding to 
the wavenumber is prov ided by the inverse of the wavenumber . 
Some other papers (e.g. Lee ct al. 1997; Abramcnko ct al. 2001) 
used exp (— i(k x x + k y y)) in the Fourier transform, in which the 
wavenumber is a factor of 2n larger than the one used in this article, 
and the spatial scale is provided by inverse of the wavenumber 
multiplied by 2tt. 
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Figure 2. Spatial power spectra of temperatures (green curves 
in the top panel), LOS velocities (blue curves in the top panel), 
and horizontal and vertical magnetic fields (purple and red curves 
in the bottom panel). The dashed and solid curves indicate power 
spectra before and after the instrument MTF and noise correction, 
respectively. The numbers indicate the power-law indices a ob- 
tained in the wavenumber range between k = 1.5 and 3.5 Mm" . 
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Figure 3. Square of the Modulation Transfer Functions (MTFs) 
of the instrument at 630 nm. The dotted curve represents a square 
of MTF of an Airy disk of a 50 cm diameter aperture. The solid 
curve is for an ideal telescope without any aberrations (i.e. the 
Strehl ratio is unity), which includes the central obscuration and 
spiders. The dashed curve is for a telescope with some aberrations 
whose Strehl ratio is about 0.8. 

there are some residual aberrations in the optical per- 
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formance. In this analysis, we use the MTF when the 
telescope has some aberration whose Strehl ratio is 0.8 
(dashed cur ve in Fig. El). R The MTF is similar to the 
one used in lDanilovic et alT(|2008| ) who studied contrast 
of granules with Hinode SP data. Though the MTF falls 
down to zero at around k = 5 Mm -1 (corresponding to 
200 km) as is shown in Fig. [31 the observed power spec- 
tra shown as dashed curves in Fig. [2] have significant 
power in the wavenumber range higher than 5 Mm -1 . 
The residual power in the high wavenumber range is 
mainly attributed to noises in the observed maps. The 
power due to the noises should be subtracted for cali- 
brating the MTF, otherwise higher power appears erro- 
neously in that range. When the data has a spatially 
uncorrelated white noise in a two-dimensional map, the 
noise spectrum Nf(k) is expected to be proportional to 
k in one-dimensional power spectra (see iRieutord et al.1 
2010). Thus the noise power Nf(k) of each parameter is 
estimated using the linear regression of the power spec- 
trum Ef{k) in the wavenumber range between k = 5 
and 6 Mm -1 . MTF-calibrated power spectra Df(k) are 
then obtained using the following equation (see lLee et al.l 
1997); 



D f (k) = (E f (k) - N } (k))/MTF{kf 



(7) 



The calibrated power spectra averaged in the quiet re- 
gions are shown as solid curves in Fig. [21 which show 
that the power in the wavenumber range higher than 
k = 0.1 Mm -1 is recovered by the calibration of the 
instrument MTF. The power in the wavenumber range 
higher than 4 Mm" 1 is acknowledged to be sensitive to 
how we subtract the noise power. Thus, the calibrated 
power spectra Df(k) are shown up to 4 Mm -1 in Fig. [21 

3. POWER SPECTRA AT THE GRANULAR AND 
SUPERGRANULAR SCALES 

The SP FOVs used to make the power spectra are a 
square of about 110 Mm, which contain several super- 
granules whose typical diameters are 20 - 30 Mm, and 
make it possible to cover both the granular and super- 
granular scales. The thermal power spectrum D t h{k) 
and the velocity power spectrum D vz (k) averaged in the 
quiet regions look very similar, as is shown in the top 
panel of Fig. [2] This can be naturally explained by the 
fact that there is a good correlation between the veloc- 
ity and the intensity due to the thermal convecti on (e.g. 
iKirk fc Livingstonlll968l iBeckers fc Parnelllll969l ) at the 
granular scale. Both the power spectra have a peak at 
k ~ 1 Mm" 1 , which corresponds to the spatial scale of 
about 1000 km and the granular scale. In the wavenum- 
ber range higher than the granular scale, i.e. in the spa- 
tial scale smaller than granules, the power spectra ap- 
pear to be obeying the power-law (i. e. oc k a ). When 
we derive the power-law indices a of the power spec- 
tra between k = 1.5 and 3.5 Mm" 1 , the indices are - 
3.3 and -3.6 for D t h(k) and D vz (k), respectively. Note 
that the power-law indices indicate that the slopes of the 

3 A wavefront error of the Hinode SOT was evaluated at G- 
band 430.5 nm with the Broadband Filter Imager (BFI) using the 
phase diversity method where on-focus and out of focus images 
were taken sequentially by moving a focusing lens (Y. Suematsu 
2012, in preparation). The MTF used in this study is derived on 
the assumption that the wavefront error derived in BFI can be 
applied to SP whose observing wavelength is 630 nm. 



power spectra at the scale smaller than granules are sig- 
nificantly steeper than k~ 5 ^ 3 which is predicted by the 
Kolmogorov's law for isotropic turbulences. The power- 
law in dices are much steeper th an the ones in iMullerl 
(1989) and lEspagnet et all (jl993) in which they reported 
that the slopes of the power spectra are close to Kol- 
mogorov's fc~ 5 / 3 power-law. Recent studies indicated 
power-law indices of about -3 based on a blue continuum 
image at 450 nm taken with the broadband f ilter imager 
(BFI) of Hinode SOT ([Rieutord et al.l [2010h and about 
-4 based on a broadband image taken through a TiO fil- 
ter at 705.7 nm with the New Solar Telescope (NST) o f 
Big Bear Solar Observatory (BBSO) ([Goode et al.l[2010h . 
which are not so different from our results. 

On the other hand, the power spectra do not show 
any enhancements in the lower wavenumber domain at 
0.03 - 0.05 Mm -1 corresponding to the supergranu- 
lar scale. This is expected because vertical flows as- 
soci ated with supergranul ation are slower than 50 m 
s _1 (jHathawav et al.l l2002f ) on the surface, and is more 
than one order of magnitude smaller than vertical flows 
of granulation, while horizon tal flows of supergranu- 
lation are 300 - 500 m s" 1 (|Simon & Leightonl Il96i 



lOrozco Suarez et al.l I2012D and are known to make a 
peak in power spectra at the supergranu lation scale 
([Hathaway et al.ll2000HRieutord et al.ll2010ft . 

The magnetic power spectra Df, z (k) and Dth{k) shown 
in the bottom panel in Fig. [2] are completely dif- 
ferent from D t h(k) and D vz (k). Both the horizontal 
and vertical powers distribute over a broad wavenum- 
ber range. There are weak enhancements at the gran- 
ular scale (k ~ 1 Mm -1 ) and at the supergranular 
scale (k ~ 0.03 - 0.05 Mm" 1 ), in both D bz (k) and 
Di,h(k), which suggests that the granular and super- 
granular scales are typical ones characterizing the mag- 
netic field distribution on the solar surface. As we men- 
tioned above, there are no significant powers at the su- 
pergranular scale in the spectra of the vertical veloci- 
ties. But it is expected that the power of the horizon- 
tal velocities are strong enough to make the magnetic 
field structures, network fie lds, at the supergranular scale 
([Orozco Suarez et al.ll2012T ) . Note that there are no clear 
enhancements of the power around the proposed meso- 
granular s cale k ~ 0.1 - 0.4 Mm" 1 , corresponding to 2.5 
- 10 Mm, ([November et al.lll981tlChou et al.lll992l ) in all 
the power spectra sh o wn in Fig. [21 which is co n sisten t 
with lHathawav et al.1 (pOOOj) and IRieutord et"aTI (|2010( ) . 
There are some residual magnetic powers in the inter- 
mediate is£ale_b^tweeji_y^grairular and supergranular 
scales. Ilshikawa k, Tsunetal ([2010D suggested that spo- 
radic appearance of granular-scale horizontal magnetic 
fields preferentially happens around boundaries of meso- 
granules (see also lYelles Chaouche et al.l 120111 ) . Such 
mesogranular-scale magnetic structures might be related 
with the residual magnetic power in the intermediate 
scale, though we need careful analysis including horizon- 
tal velocities that we cannot obtain only with the SP 
data. 

In the wavenumber range higher than the granular 
scale, the slopes of the magnetic power spectra are clearly 
less steep than the power spectra of the thermal power 
D t h{k) and the velocity power D vz {k). The power-law 
index a of Db z (k) is —1.4 in the wavenumber range be- 
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tween k — 1.5 and 3.5 Mm -1 though it is hard to say 
that the power spectrum obeys the power-law only from 
this analysis. The horizontal magnetic power spectrum 
Dbhik) has its peak at 1.3 - 2 Mm -1 corresponding to 500 
- 800 km, which is slightly smaller than the granular scale 
seen in D vz (k) and D t h(k). The spatial scale of the mag- 
netic spectra is in agreement with the observed size of 
transie nt horizontal magnetic fields (Ishi kawa fc Tsunetal 
2009b). The slope of Dbh(k) looks steeper than that of 
Dbz{k) though we have to take into account relatively 
larger noises in the horizontal magnetic power. We will 
look into this in the next section. 

4. DEPENDENCE OF THE POWER SPECTRA ON THE 
UNSIGNED MAGNETIC FLUX DENSITY 

4.1. Power spectra in small FOVs 

The power spectra shown in the previous section were 
obtained in the quiet regions, thus including both inter- 
network and network regions. In this section, we analyze 
power spectra separately in internetwork and network 
regions. To do this, each SP map is divided into 9x9 
sub-regions whose FOVs are 96 x 96 pixels, correspond- 
ing to a square of about 10 Mm, as is shown in Fig. |4j 
The FOVs of the sub-regions are not enough to cover 
the supergranular scale, but are enough to include many 
granules. The sub-regions are extracted not to use pixels 
near the edges of the SP maps. A power spectrum is 
obtained in each sub-region, using the same method de- 
scribed in §2 except the number of pixels in one direction 
N . The noise subtraction and the MTF correction are 
also applied to the power spectrum in each sub-region. 
An averaged density of unsigned flux (\B Z \) is calculated 
in each sub-region, as is shown in the right panels of 
Fig. SI The power spectra obtained in the sub-regions 
are averaged when they have similar ( \B Z |) . The interval 
of (|£?z|) for averaging the power spectra is determined 
to have 50 sub-regions in each interval. All the 30 SP 
maps are used in this analysis. Some of the sub-regions 
include sunspots in their FOVs, but such sub-regions are 
not used by checking {\B Z \). 

Fig. [5] shows the thermal, kinetic, and magnetic power 
spectra obtained in the sub-regions. The power spectra 
are the ones averaged in the interval of 5.0 Mx cm" 2 
< (\B Z \) < 5.4 Mx cm" 2 , corresponding to weak 

{\B,\) < 



internetwork regions, 13.5 Mx cm < 



14.3 Mx cm 2 , corresponding to strong internetwork 
regions, 52.2 Mx cm" 2 < (\B Z \) < 70.8 Mx cm" 2 , 
corresponding to network regions, and 229 Mx cm" 2 
< (\B Z \) < 490 Mx cm" 2 , corresponding to strong net- 
work and plage regions. The shaded area of each power 
spectrum shows variation of the power spectra (±lcr) 
in each interval. It is obvious that the magnetic power 
spectra of the horizontal component Dbh(k) have large 
variations, which is due to relatively weaker signals in 
the linear polarization. The power spectra are shown 
in the wavenumber range up to k — 3.9 Mm" 1 because 
the power in the higher wavenumber range tends to be 
affected by the noise subtraction. 

Neither thermal or kinetic power spectra, D t h(k) and 
D vz (k), change much when the average unsigned flux 
density is weaker than 20 Mx cm -2 in the internetwork 
regions (top two panels in Fig. [5]). There are peaks at 
around k ~ 1 Mm" 1 corresponding to the granular scale. 



In the wavenumber range higher than k = 1 Mm -1 , i.e. 
in the spatial scale smaller than granules, the power-law 
indices are steeper than -3 in the wavenumber range be- 
tween k = 1.8 and 3.2 Mm" 1 . The results are more or 
less consistent with the power spectra obtained in §3. 
On the other hand, the magnetic power spectra clearly 
varies depending on the average unsigned flux density 
(\B Z \). When the average flux density is very small (top 
left in Fig. [5]), the shapes of the power spectra look 
similar between Dbh(k) and Db z {k). Both the spectra 
have peaks at around k ~ 1.3 Mm -1 corresponding to 
about 800 km, which is slightly smaller than the granular 
scale and comparable with the observed size of horizon- 
tal magnetic field structures, as mentioned in th e previ- 
ous section. lOrozco Suarez S z Katsukawa (2012) exam- 
ined distribution of magnetic field inclinations in inter- 
network regions at different heliocentric angles, and con- 
cluded that the inclination distribution partly resulted 
from the granular-scale loop-like magnetic features. The 
observed peak in the magnetic power spectra shown in 
Fig. [5] probably corresponds to such loop-like features, 
which are important constituents of the internetwork 
magnetic fields. The horizontal magnetic power Dbh(k) 
is also acknowledged to be larger than the vertical one 
Df, z (k) in the internetwork regions (top tw o panels in 
Fig. Et ■ This result looks consistent with iLites et aLl 
(2008) where they reported predominance of apparent 
horizontal flux in internetwork regions. The power-law 
indices are around —1.3 for both the vertical and horizon- 
tal magnetic power spectra in the wavenumber range be- 
tween k = 1.8 and 3.2 Mm -1 . Similar power-law indices 
of the mag netic power spectra were reported in pr evi- 
ous studies (|Lee et al.lll997t lAbramenko et al.ll2001| ) al- 
though they were based on magnetic power spectra at 
the spatial scale larger than granules. Our magnetic 
power spectra show positive slopes at that scale in the 
internetwork regions. Stenflo (2012) showed that a mag- 
netic power spectrum obtained with a Hinode SP data 
set exhibited a slope that was described by a power-law 
whose index was about —3.2 for spatial scales smaller 
than 1.2 Mm. The above power- law index was derived 
from the power spectrum without applying calibration of 
the instrument MTF. In our analysis, the raw magnetic 
power spectrum before applying the MTF calibration is 
about —2.7, as is shown in the bottom panel o f Fig. 2l 
which is not so different from the one in (Stenflol [2012). 

When the average unsigned flux density is large in the 
network and plage regions (bottom panels in Fig. [5]), 
the absolute vertical magnetic power increases. The in- 
crease in magnetic power is especially noticeable in the 
lower wavenumber range (k < 1 Mm" 1 ). The horizon- 
tal magnetic power also exhibits similar behavior, and 
increases in the lower wavenumber range. The slope of 
the power Dbh(k) and Df, z (k) changes the sign from pos- 
itive to negative in the wavenumber range lower than 
k = 1 Mm" 1 , which makes the peak at the granular 
scale less pronounced in the network regions. The ther- 
mal and kinetic power spectra become affected by the 
relatively strong magnetic fields in these regions. The 
kinetic power D vz (k) tends to decrease at the granular 
scale k ~ 1 Mm" 1 , which is naturally explained by sup- 
pression of the convective motion by magnetic fields. It 
is also realized that the kinetic and thermal powers tend 
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Figure 4. Sub-regions used in the power spectrum analysis in Section 4 are shown on Hinode SP maps of the temperature fluctuation 
AT (left panels) and the vertical component of magnetic flux densities B z (middle panels). Each sub-region consists of 96 X 96 pixels 
corresponding to a square of about 10 Mm. An average unsigned flux density (|Bz|) is calculated in each sub-region as shown in the right 
panels. The top panels show maps obtained on 24 April 2007 which contain only a quiet region in the FOV. The bottom panels show maps 
obtained on 28 August 2007 which contain an active region in the FOV. We do not use the regions near the boundary of the SP maps and 
sub-regions containing a sunspot in this analysis. These regions are shown as blank areas in the (|-B Z |) maps (right panels). 



to be slightly enhanced in the wavenumber range higher 
than 2 Mm -1 . The increase in thermal and kinetic power 
at the small scale can be attributed to contributions of 
faculae structures created by small-scale magnetic fields. 
The suppression of the convection at the granular scale 
and the enhancements of the small-scale power makes 
the power spectra D t h{k) and D vz {k) less steep in the 
wavenumber range higher than k = 1 Mm . Surpris- 
ingly, all the thermal, kinetic, and magnetic power spec- 
tra show similar slopes in the wavenumber range higher 
than k = 1 Mm -1 when the average unsigned flux den- 
sity is bigger than 200 Mx cm -2 (bottom right panel in 
Fig. ED. 

4.2. Dependence of the power-law indices on the 
unsigned magnetic flux density 

Here we derive the power-law indices a of the power 
spectra obtained in the sub-regions. The power spectra 
used here are averaged ones in the interval of (|-B Z |), as 
described in the previous subsection. Thus we can get the 
power-law indices a as a function of (\B Z \). Because the 
power spectra have different slopes between k < 1 Mm -1 
(larger than the granular scale) and k > 1 Mm -1 (smaller 
than the granular scale), we derive the power- law in- 
dices a separately in these two wavenumber ranges. Here 



we use the wavenumber range between k = 0.15 and 
0.6 Mm -1 and k = 1.8 and 3.2 Mm" 1 . We do not 
say that the power spectra follow the power-law in these 
wavenumber ranges, but the indices are used to charac- 
terize the slope of the power spectra. The variation of 
the power spectra in each (|-B Z |) interval is taken into ac- 
count, and is reflected in errors of the power-law indices. 

Fig. [6] displays the power-law indices a as a function 
of the averaged unsigned flux density (|£? 2 |). Power-law 
indices of the thermal and kinetic power spectra D t h(k) 
and D vz {k) change similarly as a function of (|S Z |), while 
those of vertical and horizontal magnetic power spectra 
Dbz{k) and Dth(k) behave similarly as well within the er- 
ror bars. When the unsigned flux density is smaller than 
around 20 Mx cm -2 , the difference of the power-law in- 
dices is significant between the magnetic power spectra 
and the thermal/kinetic ones in the higher wavenumber 
range (smaller than the granular scale). The latter power 
spectra are steeper than the former ones, and the power- 
law indices differ by about 2. The magnetic and ther- 
mal/kinetic power spectra tend to have similar slopes 
when the unsigned flux density becomes large. All the 
power-law indices are around -2 when the average un- 
signed flux density is bigger than 200 Mx cm" 2 . 

In the smaller wavenumber range (larger than the gran- 
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Figure 5. Power spectra D t f l (k) (green curves), D vz (k) (blue curves), D(, z (fc) (red curves), and Di,f l (k) (purple curves) made in the 
sub-regions of the SP maps. Here shown are the spectra averaged in the intervals of 5.0 Mx cm -2 < (\B Z \) < 5.4 Mx cm~ 2 (top left), 
13.5 Mx cm" 2 < (\B e \) < 14.3 Mx cm" 2 (top right), 52.2 Mx cm" 2 < (\B Z \) < 70.8 Mx cm" 2 (bottom left), 229 Mx cm" 2 < (\B Z \) < 
490 Mx cm -2 (bottom right). The shaded area of each power spectrum represents the ±lcx range showing variation of the power spectra 
in each interval of (|B Z |). The values in the plots indicate the power-law indices a obtained between k = 0.15 and 0.6 Mm -1 and between 
k = 1.8 and 3.2 Mm" 1 . 

ular scale) , the power-law indices of the magnetic power- 
spectra decrease monotonically with the unsigned flux 
density (|-B Z |), as is shown in the upper part of Fig. [6l 
They change the sign from positive to negative when 
(|-B Z |) is 20 - 50 Mx cm~ 2 . The power-law indices of 
the thermal and kinetic power spectra do not change so 
much when (\B Z \) is smaller than around 100 Mx cm' 2 , 
and they decrease when (\B Z \) is bigger than around 
100 Mx cm~ 2 . The positive slopes become almost zero 
or negative when (|-B Z |) is around 300 Mx cm -2 , which 
means that the peaks at the granular scale disappear in 
the thermal and kinetic power spectra. When we extrap- 
olate the trend to the larger flux, the power-law indices 
in the smaller and larger wavenumber ranges might be- 
come similar. This suggests that the power spectra might 
be represented as a single slope in such large magnetic 
fluxes, though this is out of the scope of this article. 

5. SUMMARY AND DISCUSSION 

We showed the power spectral analysis of the thermal, 
kinetic, and magnetic structures on the solar surface with 
the data set taken with the Hinode SP, including calibra- 
tion of the instrumental effects. We showed in §3 that 
the kinetic and thermal power spectra had a peak at the 
granular scale while the magnetic power spectra had a 
broadly distributed power in various spatial scales, but 
had two weak peaks at both the granular and super- 
granular scales. We did not find enhancements of the 
magnetic power at the mesogranular scale though there 



were residual powers in between the granular and super- 
granular peaks. In §4 the power spectral analysis was 
done in the small sub-regions to distinguish the inter- 
network and network regions. We examined the slopes 
of the power spectra using the power-law indices, and 
compared them with the unsigned magnetic flux density 
averaged in the sub-regions. In the wavenumber range 
higher than k = 1 Mm -1 , the thermal and kinetic power 
spectra exhibited steeper slopes than the magnetic ones 
in the internetwork regions, and the power-law indices 
differed by about 2. The slopes of the thermal and ki- 
netic power spectra became less steep when the unsigned 
flux density increased, and the power-law indices of all 
the thermal, kinetic, and magnetic power spectra became 
similar. 

In the internetwork regions, where the averaged den- 
sity of the unsigned flux is very small, magnetic field 
structures are passi vely adv e cted b y the convective mo- 
tion. According to iMoffata ()196lD . when the magnetic 
field is proportional to (V • v) or (V x v) in such weak 
field regions, its power spectrum is expected to vary as 
Df,{k) cx k 2 D v (k), where and D v (k) are mag- 

netic and kinetic power spectra including all the three 
components. This relationship predicts that the slope 
of the magnetic power spectrum is less steep than the 
kinetic one, and the difference of the power-law indices 
is 2, which is consistent with the observed difference of 
the power-law indices in the higher wavenumber range 
(k > 1 Mm ) in the internetwork regions. In this anal- 
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Figure 6. Power-law indices a of the power spectra as a function 
of the averaged density of unsigned flux { | B z | ) . The power law in- 
dices are obtained in the two wavenumber ranges, between k =0.15 
and 0.6 Mm -1 (corresponding to the scale larger than granules) 
and between k =1.8 and 3.2 Mm -1 (corresponding to the scale 
smaller than granules). The green, blue, red, and purple symbols 
represent the power-lawn indices of D t f l (k), Di, z (k), Dbz(k), and 
DbhiM), respectively. The vertical error bars show ±1<t range of 
the power-law indices when we take into account the variation of 
the power spectra in each {|-B Z |) interval. 

ysis, we could not obtain the kinetic power spectrum in 
all the three components, but only in the vertical com- 
ponent. We need the horizontal velocity components 
to fully justify the scenario. When the unsigned flux 
is large, the magnetic structures are no longer passive 
to the convective motion. In this regime, magnetic and 
velocity fields strongly interact with each other, which 
possibly produces similar slopes between the magnetic 
power spectrum and the kinetic one. 

The total magnetic energy is provided by integra- 
tion of the power spectra in the wavenumber domain, 

f k 1 [D{, z (k) + Dbh(k)] dk, where fco an d fci are the lower 
and upper ends of the wavenumber range for the integra- 
tion. When both the power spectra are represented as 
a power- law (oc k a ), the total magnetic energy is com- 
pletely different in the case of between a < — 1 and 
a > —1. When a < — I (steeper than the total 

magnetic energy mainly comes from the power at the 
lower wavenumber, D bz {k ) + Dbh(k ). On the other 
hand, the total magnetic energy mainly comes from the 
power at the upper wavenumber, -Dbz(fci) + -D^fci), 



when a > — 1 (less steep than The power- law 

indices of the magnetic power spectra obtained in the in- 
ternetwork regions are close to or smaller than — 1 in the 
wavenumber range higher than k = 1 Mm . This means 
that the total magnetic energy mainly comes from cither 
the granular-scale (k ~ 1 Mm" 1 ) magnetic structures 
or both the granular-scale and smaller ones contributing 
evenly. 

Recent numerical simulations of the surface magneto- 
convection suggested that turbulent fields were created 
by small-scale dynamo action, and a significant amount 
of unsigned flux was not detected at the spatial res- 
olution of the Hinode SP because of cancellation of 
polarization signals from opposite magnetic polarities 
(Vogl er fc Schusslerl [20071: IPietarila Graham et all 120091 
120101) . These numerical simulations indicated positive 
slopes of the magnetic power spectra toward the small 
scales even at the spatial scale smaller than granules. 
The slopes do not agree with the magnetic power spec- 
tra presented in this article. If we extrapolate the ob- 
served slope fc -1 - 3 of the magnetic power spectra up 
to 50 Mm -1 , an order of magnitude higher than the 
Hinode SP resolution, the magnetic power in the ob- 
served wavenumber range between 1 and 4 Mm -1 is 
about 50 % of the power between 1 and 50 Mm -1 , 
which is not as small as the percentage suggested by 
IPietarila Graham et al.l (|2009[ ) . Because we cannot argue 
behavior of the power spectra in unresolved scales based 
on these data, it is critical to use high-resolution data 
to extend th e curre nt a nalysis toward such sma ll scales. 
iGoode et all 1)20100 and lAbramenko et al.l (|2012D showed 
kinetic power spectra of transverse velocities down to 
the 0.1 Mm scale based on high-resolution imaging data 
taken with the 1.6 meter NST of BBSO, but they did 
not have magnetic field measurements yet. One possibil- 
ity is to use magnetogram data taken with the Imaging 
Magnetograph experiment (IMaX) on board the 1 meter 
Sunr ise balloon-borne observatory ()Martmez Pillet et al.l 
I2011D to extend the power spectra down to the 0.1 Mm 
scale. iStenflol (|2012f ) argued possible existence of a flux 
tube population at the spatial scale close to or smaller 
than the resolution limit of the Hinode SP using a deep 
mode data set. Because the subtraction of the noise spec- 
trum becomes critical to derive the power in the higher 
wavenumber range, we could not argue the power at the 
spatial scale below 0.25 Mm using the normal mode data 
in this article. It is important to use polarimctric data 
with a higher signal-to-noise ratio for studying the power 
at the spatial scale as close as the instrumental resolution 
limit. 
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